Entanglement control via reservoir engineering in ultracold atomic gases 
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We study the entanglement ol two impurity qubits immersed in a Bose- Einstein condensate (BEC) 
reservoir. This open quantum system is particularly interesting because the reservoir and system 
parameters are easily controllable and the reduced dynamics is highly non-Markovian. We show 
how the model allows for interpolation between a common dephasing scenario and an independent 
dephasing scenario by simply modifying the wavelength of the superlattice superposed to the BEC, 
and how this influences the dynamical properties of the impurities. We demonstrate the existence of 
very rich entanglement dynamics correspondent to different values of reservoir parameters, including 
phenomena such as entanglement trapping, entanglement sudden death, revivals of entanglement, 
and BEC-mediated entanglement generation. In the spirit of reservoir engineering, we present the 
optimal BEC parameters for entanglement generation and trapping, showing the key role of the 
ultracold gas interactions. 
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I. INTRODUCTION 

Ultracold gases have recently emerged as an exciting 
playground for the simulation of complex many-body sys- 
tems. The great level of control and cooling over neutral 
atoms in an optical lattice has opened new avenues for the 
study of quantum magnetism, disordered systems, long 
range interactions, supersolidity, and the effect of non- 
Abelian fields, just to mention a few examples (see, for 
example, Refs. [HE])- A particularly rich area of research 
at the point between ultracold atoms, open quantum sys- 
tems and quantum information is the study of quantum 
reservoir engineering. 

Mesoscopic quantum systems can be viewed as a spe- 
cial form of environment: due to their low temperature 
and relatively small size, their quantum coherence leads 
to important effects. A Bose-Einstein condensate (BEC) 
is a typical example. Previous works have focused on the 
interaction of localised impurities immersed in a BEC 
others have studied the case of a lattice of impuri- 
ties interacting with a BEC which leads to emission of 
phonons and dissipation for the impurities @]. The col- 
lective dephasing of a two spatial states impurity was 
considered in [5], and the dynamics of a single impurity 
in a BEC was recently investigated in (6j-[8|. Moreover, 
experiments realising the dynamics of a single impurity 
immersed in a BEC [§1 HH] have recently been reported. 
All these studies and experimental progress open the way 
to applications in quantum information processing that 
were, until recently, unrealistic. 
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In general, the presence of the ultracold gas surround- 
ing the impurities gives rise to non-dissipative decoher- 
ence due to collisions between the BEC atoms and the 
impurity atoms. The literature on both entanglement 
decay and entanglement generation in the presence of 
dephasing environments has focused on certain paradig- 
matic models of environment with specific types of spec- 
tral densities [TTl - fT3] . In this paper we investigate for the 
first time a situation for which the entanglement dynam- 
ics can be modified by changing both the form of the 
spectrum and the distance between the impurities. In 
fact, in our system a modification of the BEC scattering 
length is reflected in the form of the reservoir frequency 
spectrum, leading to a change in its Ohmicity charac- 
ter . As a consequence we can investigate the optimal 
reservoir parameters to preserve both entanglement and 
discord as well as the optimal conditions for entangle- 
ment generation from a separable state. 

From a fundamental point of view our results pave 
the way to new experiments on basic aspects of open 
quantum systems. Indeed, they provide a physical exam- 
ple of an experimentally realisable bipartite open quan- 
tum system whose exact, and therefore non-Markovian, 
dynamics is described by a time-local master equation. 
This is important because, together with optical imple- 
mentations [H], these systems could be used as testbeds 
for experimental verification of fundamental theorems of 
open quantum systems, thus giving an essential contri- 
bution to the thriving field of non-Markovian quantum 
dynamics [IB]. More specifically, the two- impurities in 
BEC we discuss in this paper could be used to verify the 
simplest nontrivial generalisation of the Lindblad-Gorini- 
Kossakowski-Sudarshan (LGKS) theorem for commuta- 
tive dynamical generators, as we briefly discuss in Sec. 
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II. 

The structure of the paper is the following. In Sec. II 
we recall the microscopic derivation of the effective pure 
dephasing model, we derive the master equation for the 
two impurities, and discuss its interest for fundamental 
tests of open quantum systems dynamics. In Sec. Ill we 
discuss the dynamics of both entanglement and quantum 
discord for different reservoir parameters. In Sec. IV we 
show the onset of entanglement generation and its de- 
pendence on the distance between the qubits and on the 
scattering length. Finally, in Sec. V we present conclu- 
sions. 



following Hamiltonian [5] 




FIG. 1. (Color online) Schematic view of the proposed setup. 
Two impurities are each separately trapped in a deep dou- 
ble well potential, such that tunnelling between wells is sup- 
pressed. In this regime each atom occupies the two states 
localized in each well. The two impurities are immersed in 
a three-dimensional BEC with which they interact via cold 
collisions. We denote by L and D the separation between the 
two wells in an impurity and the separation between the two 
double-wells, respectively. 



II. THE MODEL 



The system we consider consists of two impurity 
atoms immersed in a Bose-Einstein condensate reservoir 
trapped by a large harmonic potential. In contrast to 
similar proposals in which atoms were confined in a tight 
harmonic trap [3JH], in the present model, sketched in 
Fig. [lj each atom is trapped in a double well potential. 
The two minima are separated by a distance L, with the 
double wells separated by distance D. Our qubit consists 
of the presence of the atom in the left or right well of the 
double well. The combined system is governed by the 



H = H A + H E 
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is the Hamiltonian of the atom impurity with mass tua, 
where ^(x) is the corresponding atomic field opera- 
tor [22], Va(x) is the double well potential formed by 
the optical lattice, 



Hr 



dxd^x) 



A + Vb(x) + ?f $t( x )$( x )l $( x ) 
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is the Hamiltonian for the BEC atoms of mass m B , 
is the corresponding field operator satisfying bosonic 
commutation relations, Vb(x) is the BEC trapping po- 
tential which we assume to be constant and gs — 
Anh 2 aB/mB is the boson-boson coupling constant, with 
as the BEC scattering length. We denote by 



Hab — 9ab J dx*t( x )$t( x )$( x )$( x ) 



(4) 



the interaction Hamiltonian of the impurity atoms with 
the condensate, and gAB = 2-kTi} 'a ab / ab is the cou- 
pling between the impurity atoms and the condensate 
gas. We defined the reduced mass mAB = rn A mB / (m^l 
ms) and the scattering length ciab corresponding to 
impurity-boson collisions. We now make a number of 
approximations leading to a simplified model. First we 
assume deep double wells, meaning that the impurity 
atoms can only be found in the lowest Wannier states, 
\L) and \R), of the double wells corresponding to the 
atom localised in the left and right well respectively and 
we neglect tunnelling across the barrier separating the 
two wells. This allows one to rewrite Hamiltonian (Til) 
using the pseudo-spin operator a z = \L) (L\ — \R) (R\. 
The second approximation valid for a weakly interacting 
Bose gas, consists of expanding the field operator as 



<D(x) 



AT $o(x) + [ M k(x)c k - i£(x)c\ 

k 



(5) 



where $o( x ) is the classical condensate wave function and 
the operators Ck correspond to the Bogoliubov modes of 
the BEC (see for example [33] ). Using this assumption, 
and following the derivation of Ref. [5] we arrive to an 
independent boson model with an effective Hamiltonian: 



h = H Bog + \ E - 



,2=1,2 



i=l,2 



^L,k) C k + h.C. 



(6) 
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where 



(7) 



is the Bogoliubov Hamiltonian and Ek is the correspond- 
ing energy. The constants f2^. k and fl l L k are the coupling 
frequencies of the left and right localised states of qubit i 
respectively to the k-th Bogoliubov mode, and h.c is the 
Hermitian conjugate. The form of the constants Q de- 
pends on the overlap integrals of the impurities Wannier 
functions with the wave functions associated to the Bo- 
goliubov modes. Their complete form and their deriva- 
tion can be found in [5]. We assume a constant BEC 
wave function <J> that is not modified by the presence of 
the impurities. This means that the terms of the Hamil- 
tonian responsible for the Gross-Pitaevskii Equation do 
not enter into Eq. [6j Finally a\ is the Pauli operator 
of qubit i. Further details about possible experimental 
realisations are given in Appendix |A"} 

We assume an initially factorized state of the environ- 
ment and the two qubits: 



g(0) = p(0) (8 t b 



(8) 



where p(0) is the initial (possibly correlated) density ma- 
trix of the two qubits that we will specify later and 
tb = cxp(—Hb/kbT)/Z is the condensate excitations 
density matrix assumed to be in equilibrium at tem- 
perature T (fig is the Boltzmann constant and Z = 
Trexp(-H B /K B T)). 

Since we are interested in the two-qubit dynamics only, 
we trace out the state of the environment and compute 
the reduced density matrix after time t has passed: 



p(t) = Tr B g(t). 



(9) 



The populations of the density matrix do not change 
with time. This is a consequence of the form of the 
impurity-condensate density-density interaction which 
induces pure dephasing in the impurity system. On 
the other hand the coherences, represented by the off- 
diagonal terms, are governed by a non trivial dynamics: 



-r w (i)+in y (t) 



(10) 



The quantities (f ) and (t) are real functions of time 
and describe decay of coherences, that is, proper decoher- 
ence, and phase shifts, respectively. The precise form of 
these functions was computed in [5] and here we only 
report the form for the T factors. Here we focus on a 
three-dimensional environment . 

The off diagonal elements each decay according to one 
of three decoherence parameters. The first: 
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where f3 = l/(regT), is responsible for the decoherence of 
each qubit individually. Therefore it appears in matrix 
elements such as P{ll},{lr} or P{ RL },{rr}- 

The other two factors entering in matrix elements such 
as P{ll},{rr\ or P{lr},{rl}, in which the collective co- 
herence of the two qubits is involved, are 



T ± (t)= 2 -^J Q 00 dkk 
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(12) 



where 8(t) represents the cross talk between the qubits 
via the reservoir. For very short times, due to their spa- 
tial separation, each qubit decoheres independently as 
shown in [SJ[S]. For larger times, determined by the ratio 
of the distance over the typical speed of sound, the three 
parameters become distinct and the effects of the com- 
mon reservoir become apparent. Increasing the distance 
between the two impurities allows for simulation the dy- 
namics of two qubits in local dephasing environments. 
Hence, by changing D, that is, the wavelength of the su- 
perlattice, one can pass from a common environment to 
an independent environment scenario. 

Starting from Eq. ( 10 ) a straightforward calculation 



allows us to prove that the exact master equation for the 
system can be written in the following time local form, 



d P (t) 1 • 



(a a z - a b z )p(a a z - a\) - -{{a a z - a b z )(a a z - o\\ p} 



(13) 



where the time-dependent dephasing coefficients are the may attain negative values, indicating that the dynamics 

time derivatives of the decoherence factors given by is non-Markovian or, more precisely, that the dynami- 

Eqs. (12 1, and where the jump operators are given by cal map is non-divisible. It is worth stressing that the 

cr° ± a° z . One can verify that the dephasing coefficients general conditions under which a non-Markovian master 
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equation can be written in time-local form are not yet 
known. To the best of our knowledge this is the first 
example of a two-qubit master equation derived from a 
microscopic model which is proven to be in time local 
form. 

Our system is a simple example for which the gener- 
ators of the dynamics form a commutative algebra. For 
these systems one can generalise the Lindblad-Gorini- 
Kossakowski-Sudarshan theorem [IT] proving that the 
dynamics is completely positive if and only if the de- 



T = 10"°K 



T= 10"'K 



coherence factors of Eqs. (12) are positive at all times 



[18) . The physical implementation and verification of 
our model would allow then to verify the generalization 
of the LGKS theorem for the simplest nontrivial com- 
mutative open quantum system. Both the single qubit 
dephasing case and the dephasing of two qubits in lo- 
cal environments, indeed, have been shown to posses a 
random unitary representation |19j . 
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FIG. 2. (Color online) Entanglement phase diagram for p^ as 
a function of scattering length, as in units of am,, and initial 
state parameter, c, for D = 5A and a) T = 10 _S K and b) T = 
10 -7 K. We distinguish three phases: entanglement trapping 
(dark blue, top right), entanglement sudden death (white, 
bottom left) and entanglement revivals (medium blues) (see 
main text). 



III. CORRELATIONS FOR INITIALLY 
ENTANGLED STATES 

Let us consider the dynamics of entanglement between 
the impurities, and its dependence on the system and 
reservoir externally controllable parameters. Due to their 
computational simplicity combined with their potential 
for entanglement and non-classical behaviour, we con- 
sider as initial states Werner states of the form 



p+ =c|$+) 



1 - c. 



ce[o,i], (14) 



and 



Pw = c\*+)^ + \+ l -^\ ce[0,l], (15) 

where |$+) = (\LL) + \RR))/y/2 and |*+) = (\LR) + 
\RL))/y/2 are maximally entangled states and I is the 
4-dimensional identity matrix. As an example, in Ap- 
pendix [S] we discuss how to prepare the two qubits in 
the state given by Eq. (14 1. As for the generation of 
the second class of initial states here considered, given 



by Eq. (15), we note that all Bell states are equivalent 
up to local operations, i.e., if we are able to create one 
then we can generate all others, therefore the procedure 
sketched in Appendix [A] can be generalised accordingly. 
Werner states have applications in quantum information 
and quantum teleportation [53J [S3] . Additionally, they 
interpolate between maximally entangled and separable 
states, that is, they can have any value of concurrence 
between and 1, and therefore are ideal states for in- 
vestigating entanglement dynamics in our system. For 
c < 1/3 the state is separable; above that, the state ini- 
tially has non-zero concurrence which will time evolve 
as 



(*) 



0, ce 



-r T («) _ 



1 
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FIG. 3. (Color online) Stationary concurrence as a function of 
relative scattering length, as/aRb, for D — A/2, T = 10 _8 K 
and for initial states p^ (green dashed line) and p^y (blue 
solid line). 



This formula shows that, if the system is initially pre- 
pared in a Werner state, the interaction with the con- 
densates cannot increase the initial entanglement and in 
general leads to a loss of concurrence. However, as we will 
see in the following, by manipulating the system-reservoir 
parameters we can control entanglement dynamics and 
maximise the amount of stationary entanglement in the 
system. From Eq. (TTsJ) we also see that, for an initial 
state of the form pyy (p^y) the only decay rate which 
enters the dynamics is the collective decay rate r_(t) 

(r+(t)). 

It is worth noticing that, for the states here consid- 
ered, Cf^(i) does not depend on the phase shifts Hy(t). 
The behaviour of the concurrence depends primarily on 
three factors: the initial state parameter, c, the scattering 
length of the condensate, ag, and the distance between 
the qubits, D. By varying these parameters, it is pos- 
sible to obtain a wide range of entanglement dynamics 
for long times. Specifically three different types of en- 
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tanglement dynamics can be observed: (i) sudden death 
of entanglement, (ii) sudden death followed by revivals, 
(iii) entanglement trapping. 

Fig. [2] (a) shows the distribution of the types of entan- 
glement behaviours for two qubits at a distance D = 5A, 
where A is the wavelength of the optical lattice, for vary- 
ing initial parameters, c, and scattering lengths, as, mea- 
sured in units of the natural 87 Rb scattering length a Kb, 
for the initial state of Eq. ( fl4"| . The white region corre- 
sponds to entanglement sudden death, the medium blue 
region corresponds to sudden death followed by revivals 
and the dark blue region corresponds to entanglement 
trapping. It is worth stressing the extreme sensitiv- 
ity of the dynamics to the initial state, in the interval 
0.38 < c < 0.44. We recall that the state is separable 
for c > 1/3 ~ 0.33. Entanglement sudden death takes 
place for small values of initial entanglement, that is, for 
0.33 < c < 0.38, while if the initial entanglement is suffi- 
ciently high the phenomenon of entanglement trapping, 
a nonzero value of stationary entanglement, will occur. 
This is the case for any value of c > 0.425, independently 
on the value of the scattering length, hence we can con- 
clude that in this system entanglement trapping is dom- 
inant. A similar entanglement phase-diagram exist for 
initial states of the form of Eq. ( p~5[ ) . 

Entanglement trapping originates from the fact the the 
time dependent dephasing rates appearing in the master 
equation (13) go to zero after a finite time. Hence de- 



phasing stops and so does entanglement loss. This is a 
strongly non-Markovian phenomenon and never occurs 
for systems described by master equation with positive 
constant decay rates. 

In between the sudden death and entanglement trap- 
ping regions there are two small regions where entan- 
glement exhibits periodic revivals which may or may 
not result in residual entanglement. It is worth notic- 
ing that, even when temperature is increased, there is 
still entanglement trapping for the majority of cases, as 
shown in Fig. [2] (b). The effect of the temperature is 
to enhance entanglement decay and enlarge the param- 
eter space for which entanglement sudden death occurs. 
In fact, for typical experimental temperatures of the or- 
der of T = 10~ 8 K to T = 10 _7 K entanglement trapping 
may only take place for sufficiently strongly interacting 
gases. In general, an increase in the scattering length of 
the ultracold reservoir will favor entanglement trapping, 
allowing for this phenomenon to occur for a wider class 
of initial states. Moreover, as we can see from Fig. [3] the 
stationary entanglement increases for increasing values 
of as, for both classes of initial states. This conclusion 
holds for any value of the distance D between the impu- 
rities. Hence we conclude that more strongly interacting 
gases are optimal for reaching higher values of entangle- 
ment trapping. 

We now investigate the change in stationary entangle- 
ment when increasing the qubits separation, for a fixed 
value of the scattering length. Figure [4] shows and 
as a function of D/L. Initially, there is an increase 
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FIG. 4. (Color online) Steady state value of concurrence, 
C^ nal , as a function of distance between qubits for initial 
states (blue solid) and (green dashed) for as = am- 



or decrease, respectively, of the entanglement, settling to- 
wards a steady value as the qubits get further apart. The 
different initial behaviour of stationary entanglement for 
the two classes of initial states can be explained as fol- 
lows. When immersed in a common environment the ini- 
tial states pyy, which are mixtures of the subradiant state 
|$ + ), are sub-decoherent [5J. In our model an increase 
in the distances corresponds to a vanishing effect of the 
cross talk term, 5, and hence to a transition to a model 
of local dephasing for the impurities. For locally dephas- 
ing impurities sub-decoherence does not occur. This ex- 
plains why an increase in the distance causes a decrease 
in the trapped entanglement for initially sub-decoherent 
states. The opposite holds for the other class of initial 
states p^, which are super-decoherent in a common en- 
vironment. In general, for large D, the cross talk term, 
5, vanishes and the qubits behave as single qubits in in- 
dependent reservoirs, eliminating the difference caused 
by their initial states. In terms of reservoir engineering, 
states p^r are best protected from decoherence when im- 
mersed in a common environment, that is, for shorter 
values of the impurity separation due to smaller super- 
lattice wavelength, while states p^ maximally retain the 
initial entanglement for longer impurity separations, that 
is, larger superlattice wavelengths. 

Until now, we have analysed the reservoir parame- 
ters maximising the value of residual entanglement. We 
now look at the changes in the entanglement dynam- 
ics for different initial states and impurity distances. 
We focus our attention on the most interesting region 
of system-reservoir parameter space where the entangle- 
ment is most sensitive to the initial conditions. Fig. [5] (a) 
shows the concurrence for two qubits at D = 2L = A/2, 
that is, in adjacent pairs of sites of an optical lattice, for 
three values of the initial state parameter, c. As c in- 
creases, we can see the concurrence move from zero value 
to a finite value in the steady state. Fig. [5] (b) shows the 
same for qubits located at D — 5A. At these distances the 
cross talk term, S, picks up oscillations and the resulting 
concurrence shows more varied behaviour. In addition, 
the same initially entangled state will result in different 
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FIG. 5. (Color online) Concurrence, Ci, for a) D — A/2 for 
initial state parameters c = 0.34 (blue solid), 0.341 (green 
dashed) and 0.342 (red dotted), and b) D = 5A for initial 
state parameters c = 0.339 (purple dotted), 0.3395 (orange 
dashed) and 0.34 (blue solid). 
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FIG. 6. (Color online) Concurrence for the initial Werner 
state (blue solid) and discord (green dashed) as a function 
of time for as = a-m and D — 10L for initial parameter 
c = 0.3395. 



concurrences for different values of D. In particular, for 
the same initial state, c = 0.34, the steady state value of 
the concurrence changes from zero in Fig. [5] (a) to a small 
but non-zero value for Fig. [5] (b) due to the increased dis- 
tance between qubits. 

For the sake of completeness, we conclude our analy- 
sis of correlations dynamics by looking at the quantum 
discord of the impurities, as an example of a quantum 
correlation that exceed those captured by entanglement 
|21j . We consider the original definition of discord: 



Q(p)=l(p)-C(p), 



(17) 



where I(p) is the quantum mutual information and C(p) 
quantifies the classical correlations. As such, this quan- 
tity describes correlations beyond those ascribable to 
classical physics and highlights the presence of quantum 
effects even in cases where there is no entanglement. For 
Werner states of the form outlined above, the discord 
can be calculated using a simple algorithm |26i . In Fig. [6] 
we compare the discord and concurrence as a function of 
time. Here the concurrence shows an initial loss of entan- 
glement, followed by periods of revival. The behaviour 
of the discord follows the same general trends as the con- 
currence, having peaks for the same instants of time. For 
this class of states, this is probably due to the simple de- 
pendance of both quantities on qubits correlations {(J z ct z ) 
and (cr x a x ). An important difference, however, is that 
in the time intervals in which entanglement temporar- 
ily disappears, there are still oscillations present in the 
discord showing that the state still displays non-classical 
behaviours. This is not surprising as, in general, it is 
known that quantum discord is more robust against de- 
coherence than entanglement and, for example, does not 
exhibit the phenomenon of entanglement sudden death. 



IV. INITIAL PRODUCT STATES AND 
ENTANGLEMENT GENERATION 

In the previous section, we have considered the effect 
of the environment on states that are initially entangled 



and have shown that the presence of a carefully chosen 
shared environment protects the system from loss of en- 
tanglement. However, one key advantage of a shared en- 
vironment is its ability to generate correlations between 
initially separable states. In analogy to (T3], here we con- 
sider an initial product state, 

|V(0)> = i (|LL) + \LR) + \RL) + \RR)) , (18) 

and evolve it in time. In contrast to Werner states in 
which only one of the decay rates, T±(t), would appear, 
for this initial state the dynamics is dictated by all three 
decay rates, including r . Moreover in this case the en- 
tanglement crucially depends on the phase shifts Hy (i). 
In order to measure any generated non-classical correla- 
tions developed in the system, we calculate the concur- 
rence of the resulting state. Fig. [7] shows the concur- 
rence as a function of time for this initial state. The 
concurrence shows a periodic oscillation between its zero 
value and a value close to its maximum 1, for temper- 
atures of T = 10 _8 K. If the temperature is increased 
to T = 10 _6 K, the periodic generation of entanglement 
remains, but the maximum attainable concurrence is re- 
duced. 

The origin of this entanglement and the specific choice 
of the particular initial state ( 18 1 can be explained by 



the fact that the BEC induces an effective interaction 
Hamiltonian proportional to a z a z . This interaction cre- 
ates conditional shifts between the two qubit impurities 
which, for a state of the form (18), correspond to an 



entangling operation. Therefore for other product states 
such as | Li) or \LR) we do not observe any entanglement 
generation at all. 

We now explore how the generated entanglement can 
be maximised by varying reservoir parameters. Fig. [8] 
shows the dependence of both the maximal entangle- 
ment and the generation time when one changes either 
the scattering length as or the distance D between the 
impurities. As the entanglement generation crucially re- 
lies on the phase factors, and as these quantities have 
negligible effect for increasing distances D, one sees that 
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the maximum concurrence decreases with increasing dis- 
tance while the generation time increases. On the other 
hand, the plots show that a more strongly interacting 
BEC leads to higher values of entanglement, that is, once 
more the optimal reservoir configuration is for higher val- 
ues of as ■ The price to pay is a greater value of the gen- 
eration time, indicating that one has to wait longer for 
the system to attain the highly entangled state. 

Apart from fundamental implications for quantum 
reservoir engineering, our setup may have applications 
in optical lattices for the creation of long distance en- 
tanglement between impurities trapped in a superlattice. 
One does not need to manipulate the form of the im- 
purities' potential by, for example, lowering and raising 
the barrier in a double well. One needs to control the 
time interval that impurities and BEC interact by mov- 
ing the impurities away from the BEC or by lowering to 
zero the interaction strength qab by means of Feshbach 
resonances. 
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FIG. 7. (Color online) Concurrence for the initial state of 
Eq. ( 18 \ as a function of time for T = 10 _8 K (solid blue) and 
T = 10" 6 K (green dashed). 



V. CONCLUSIONS 

Atoms immersed in Bose-Einstein condensates pro- 
vide an ideal system for investigating and probing many- 
body dynamics thanks to their interaction with the quan- 
tum excitations of the condensate. We have shown that 
the combination of the double well qubits and a Bosc- 
Einstein condensate reservoir offers several advantages 
for quantum information processing and storage. For ini- 
tially entangled Werner states, there is an initial loss of 
entanglement to the environment, but the presence of the 
common environment reduces, and in some cases partly 
reverses, the loss of correlations from the system to the 
environment. Moreover, for certain initially separable 
states, the BEC-induced effective interaction can gener- 
ate entanglement between the two distant qubits. In the 
spirit of reservoir engineering, we have demonstrated how 
we can manipulate both the entanglement dynamics and 
the residual entanglement by modifying system-reservoir 
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FIG. 8. (Color online) a) Cw ax , the maximum concurrence 
generated for the initial state of Eq. (18 1, as a function of 
scattering length as I am for D = A/2 and T = 10~ 8 K. b) 
tmax , the time at which the maximum concurrence of (a) oc- 
curs, 
c) 



state of Eq. ( 18 1, as a function of distance between qubits, D 
for as = am and T = 10 _8 K. d) t max , 
maximum concurrence of (c) occurs. 



the time at which the 



parameters and we have shown that more strongly in- 
teracting BECs are optimal for both entanglement gen- 
eration and entanglement trapping. Moreover, we have 
shown that the dependence on the distance between im- 
purities is non trivial. If one aims at prevention of 
entanglement, then different classes of initial Werner 
states have opposite requirements (D/L ~ 1 for and 
D/L 1 for p^), while entanglement generation is al- 
ways optimal for small distances. 

Finally we have shown that our results are robust 
to the effects of experimentally realistic temperatures, 
meaning that our proposal may be tested in present-day 
experimental setups. Our results pave the way to im- 
plementations of quantum communication protocols in 
arrays of quantum impurities and underline at the same 
time the potential of these systems as testbeds for fun- 
damental studies on open quantum systems in the non- 
Markovian regime. 
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Appendix A: Experimental realization 

The ideas and schemes presented in this work can be 
implemented using techniques realised in recent experi- 
ments of impurity-BEC dynamics. We consider an im- 
plementation with 133 Cs impurity atoms immersed in a 
87 Rb condensate, but our results can be extended to 



other species or a single specie BEC but with two dif- 
ferent internal states. The impurities can be trapped 
by an optical superlattice with a double well elementary 
cell. In the simplest case D = 2L = A/2. For the case 
described in Fig. [4] in which L and D are not in a simple 
ratio, one needs to use a different trapping mechanism 
for the impurities, for example arrays of microtraps [27] 
or of optical tweezers [2"B"j . 



The initial Werner states considered in Eq. (14) can 
be created using an entangling operation as the "square 
root of swap" operation realized in |29j for two internal 
states |0) and |1) of the two atoms. The square root of 
swap gate is defined as: 



VSWAP 



/ 1 \ 

LA/2 LA/2 
-1/V2 1/V2 

\ 01/ 



(Al) 



The protocol to create the Werner state is as follows: 

i) prepare the two atoms in two distinct single wells; with 
probability c create the state 1 10} and with probability 
1 — c create the maximally mixed state of the two qubits 
1/4 by applying a fast random phase to an equal super- 
position of the two states |0) and |1) for the two atoms; 

ii) apply the v / SWAP gate by bringing the two atoms 
together, let them interact and separate them again; iii) 
finally transfer the internal state of each atom into the 
left and right states of the double well, that is, |0) — > \L) 
and |1) — > \R) using a spin dependent potential and at 
the same time transforming each single well into a double 
well by raising the central barrier. To create the state 
use 1 01) in stage i). A similar procedure, not involving 



the VSWAP can be used to create state Eq. (pL8[) . 



Finally the readout can be achieved by transferring 
back the positional states \L) and \R) into the internal 
states |0) and |1) and then doing a full tomography of 
the two-qubit state. 



[1] I. Bloch, J. Dalibard and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 
[2] M. Lewenstein, A. Sanpera and V. Ahufinger: Ultracold 

Atoms in Optical Lattices: Simulating quantum many- 
body systems, Oxford University Press (Oxford, 2012). 
[3] A. Recati, P. O. Fedichev, W. Zwerger, J. von Delft, and 

P. Zoller, Phys. Rev. Lett. 94,040404 (2005); A. Klein 

and M. Fleischhauer, Phys. Rev. A 71, 033605 (2005); 
[4] A. Klein, M. Bruderer, S.R. Clark and D. Jaksch, New 

J. Phys. 9, 411 (2007) 
[5] M. A. Cirone, G De Chiara, G. M. Palma and A. Recati, 

New J. Phys. 11, 103055 (2009). 
[6] P. Haikka, S. McEndoo, G. De Chiara, G. M. Palma, and 

S. Maniscalco, Phys. Rev. A 84, 031602 (2011). 
[7] F. Mulansky, J. Mumford, and D. H. J. O'Dell, Phys. 

Rev. A 84, 063602 (2011). 
[8] S. Peotta, D. Rossini, M. Polini, F. Minardi, R. Fazio, 

arxiv: 1206.3984. 
[9] S. Will, T. Best, S. Braun, U. Schneider, I. Bloch, Phys. 

Rev. Lett. 106, 115305 (2011). 
[10] N. Spethmann, F. Kindermann, S. John, C. Weber, 

D. Meschede and A. Widera, Appl Phys B (2012) 

106:513519; N. Spethmann, F. Kindermann, S. John, C. 

Weber, D. Meschede and A. Widera, arXiv:1204.6051. 
[11] G. M. Palma, K.-A. Suominen and A. K. Ekert, Proc. R. 

Soc. Lond. A 452, 567 (1996). 
[12] U. Weiss, Quantum Dissipative Systems (1992), World 

Scientific. 

[13] D. Braun, Phys. Rev. Lett. 89, 277901 (2002). 

[14] B.-H. Liu, L. Li, Y.-F. Huang, C.-F. Li, G.-C. Guo, E.- 

M. Laine, H.-P. Breuer, and J. Piilo, Nature Phys. 7, 931 

(2011). 



[15] A. J. Leggett et al., Rev. Mod. Phys. 59, 1 (1987). 

[16] See, for example: H-P. Breuer, E.-M. Laine, and J. Piilo, 
Phys. Rev. Lett. 103, 210401 (2009); E.-M. Laine, J. 
Piilo and H-P. Breuer, Phys. Rev. A 81, 062115 (2010); 
A. Rivas, S. Huelga and M. Plenio, Phys. Rev. Lett. 105, 
050403 (2010), D. Chruscihski, A. Kossakowski, A. Rivas, 
Phys. Rev. A 83, 052128 (2011), X.-M. Lu, X. Wang and 
C. P. Sun, Phys. Rev. A 82, 042103 (2010). 

[17] G. Lindblad, Comm. Math. Phys. 48, 119 (1976); V. 
Gorini, A. Kossakowski, E.C.G. Sudarshan, J. Math. 
Phys. 17, 821 (1976). 

[18] D. Chru?ci?ski and A. Kossakowski, arXiv:1201.5987. 

[19] A. Pernice, J. Helm, and W. T. Strunz, arXiv: 1202.4280. 

[20] W. K. Wootters, Phys. Rev. Lett. 80, 2245, (1998). 

[21] H. Ollivier and W. H. Zurek, , Phys. Rev. Lett. 88, 
017901 (2001); L. Henderson and V. Vedral, J. Phys. A 
34, 6899 (2001). 

[22] The bosonic or fermionic nature of the impurity and thus 
the commutation relations fulfilled by VP are irrelevant for 
our model as there is only one impurity per double-well. 

[23] L. Pitaevskii and S. Stringari, Bose-Einstein Condensa- 
tion (Oxford University Press, Oxford, 2003). 

[24] M. Horodecki, P. Horodecki, R. Horodecki, Phys. Rev. A 
60, 1888 (1999). 

[25] J. Lee and M.S. Kim, Phys. Rev. Lett. 84, 4236 (2000). 

[26] M. Ali, A. R. P. Rau, and G. Alber, Phys. Rev. A 81, 
042105 (2010). 

[27] A. Lengwenus, J. Kruse, M. Schlosser, S. Tichelmann, G. 

Birkl, Phys. Rev. Lett. 105, 170502 (2010). 
[28] C. Muldoon, L. Brandt, J. Dong, D. Stuart, E. Brainis, 

M. Himsworth, A. Kuhn, arXiv:1109.0657 
[29] M. Anderlini, P. J. Lee, B. L. Brown, J. Sebby-Strabley, 

W. D. Phillips, J. V. Porto, Nature 448, 452 (2007). 



